Manganites at Low Temperatures and Light Doping: Band Approach and Percolation 
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A tight-band model is employed for the e2 S -orbitals in 
manganites. It is shown that a large intraatomic Hund's cou- 
pling, Jh and the resulting double-exchange mechanism lead 
to antiferromagnetic ordering along one of the cubic axis sta- 
bilized by the cooperative JT effect which further decreases 
the band energy of electrons. As a result, LaMn03 is a band 
insulator built of 2D ferromagnetic layers. The critical con- 
centration x c ~ 0.16, for onset of ferromagnetic and metal- 
lic behavior at low temperatures in Lai-^Sr^MnOa and the 
phase transition are treated in the percolation approach. 

PACS numbers: 75.10Lp; 75.30.-m; 75.50-y 

Current efforts in studies of the "colossal magnetore- 
sistance" (CMR) in the manganites, Ri_ x B x Mn03 (usu- 
ally, R=La, B=Ca, Sr; see e.g. [1-3] for a review and 
references) are, naturally, focused on phenomena occur- 
ring in the vicinity of the metal-insulator transition tem- 
perature, T*, which for x = 0.2 — 0.3, reaches the room 
temperature range. There are no doubts that both the 
Jahn- Teller (JT) distortions and Zener's double exchange 
(DE) mechanism [4,5] are two key ingredients of CMR. 

The transition between the paramagnetic and conduct- 
ing ferromagnetic phases is treated in [6] as a localization- 
delocalization "crossover" driven by thermal fluctuations 
of the local JT-modes. In this paper we focus on the 
analysis of phenomena occurring at low T. In this region 
the approach should be different. Indeed, the two end 
members of the series, LaMnOa and CaMn03,are both 
antiferromagnetic (AF) insulators. "Doping" of LaMn03 
by divalent atoms results in metallic conductivity only 
above some critical concentration, x c ~ 0.16. 

Low temperature phases of manganites are habitually 
interpreted in terms of localized orbitals [7, 8], more often 
that not in the modified Hubbard model (see. e.g. [9]). 
Although it is widely accepted that electron-electron cor- 
relations play the key role in low temperature behavior of 
manganites, we show that their major properties can be 
understood within the framework of the band approach, 
as well, or in percolation terms. 

The following facts have to be explained: 1) the in- 
sulating state of the parent LaMnOa; 2) the A-type 
of antiferromagnetic ordering (alternating ferromagnetic 
planes); 3) the small value of the Neel temperature, TV ~ 
150K in the low x-range, while the structural changes 
occur at about 900K; 4) the localization of "holes" in- 
troduced into the Mn 3+ -subsystem at x < x c ; 5) the 
threshold concentration, x c ~ 0.16, above which metallic 



conductivity sets in at low T. 

In the band model the properties from 1) to 4) are 
brought about by the DE- and JT-mechanisms. We in- 
terpret the value x c ~ 0.16 in terms of percolation theory: 
both the cluster approach and the phase separation pic- 
ture, as suggested first for the cuprates [10], seem to give 
the same criterion [11]. 

Let us start from the A-type of the ground state of 
LaMnOa. According to our picture, appearance of such a 
magnetic state is not caused by any exchange interaction 
between localized spins on different sites. The mutual 
arrangement of the distorted octahedra and of the core 
3 ^2g-spins, Si (S=3/2), is stabilized only by gains in the 
kinetic band energy of the "conduction" e2 S -electrons. 

The single electron Hamiltonian is then of the form [6] : 
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JhS, -a + gr, ■ Q t + J e ;Q?) (1) 



Here ta+s is the nearest neighbor tunneling matrix, de- 
fined on some basis of the two-dimensional cubic repre- 
sentation, e2 g ', Qi are the active local JT-modes with the 
matrix elements on that electronic basis in terms of a 
"pseudospin" matrix, r (see e.g., in [7]; Qi arc defined 
as dimensionless parameters). The static JT deforma- 
tions and the staggered magnetization, (Si) = (— 1)'(5) 
of the A-type phase are treated below in the mean field 
approximation. 

To avoid cumbersome expressions, we first discuss the 
competition between ferro- and antiferromagnetic order 
for the one band model: 



e(p) = tcos(ap z ) +t(pj 



(2) 



Assume first that the local ( 3 i2g) spins are ordered fer- 
romagnetically. At large Jjj t the low energy band is 
shifted down, by —Ju(S). The main contribution to the 
total energy is Eh = —Jn(S)n (n is the number of elec- 
trons per unit cell) , while the kinetic energy contribution 
due to the polarized electrons is linear in t, i. 

Consider the same problem for AF-ordering along the 
z-axis (with the period 2a). It is convenient to discuss a 
more general case: 



Si = (±(S Z ),M X ) ■ (S 2 Z + Ml = S 2 ) 



(3) 



i.e. the canted AF-structure [8] for the 3 t 2g - spins. Solv- 
ing the periodic electron spectrum problem in a manner 
similar to the solution for the two local sites [5], one ob- 
tains the following four branches: 



1 



ej(p) - i(p_i_) = 
±^J 2 H S 2 +t 2 {p z ) ± 2J H | t( Pz ) \\M X \ (4) 

At Jh ^> t, t, the low energy spectrum reduces to: 

e(p)i, 2 - t(p±) - - Jh^± I i(p*) | (| M x | /S) (5) 

The second term restores the familiar expression for 
transport [5]: M x = 2S | cos^/2 | (9 is the angle between 
the adjacent local spins). 

For the AF-case (M x = 0), the Brillouin zone is re- 
duced by half, but there are now two branches in (5). As 
the result, the main term, — JhS, does not change. For 
a single band, the antifcrromagnetic order may even be 
energetically favorable, although the energy gain would 
be small, of the order of (see (4)): 

t 2 /J H < t (6) 

For two bands there are terms linear in t. The JT-band 
splitting is necessary if we wish to decrease the ground 
state energy below that one for the ferromagnetic state. 

We obtained the "bare" band spectrum of the cubic 
material with the e2 9 -electrons. For calculating tik, and 
hence, the electronic spectrum, e(p), use of the normal- 
ized basis functions of the form: 

ip! oc z 2 + ex 2 + e 2 y 2 , V>2 = i>* (7) 

(e = cxp(27ri/3)), proves to be more convenient to ac- 
count for the cubic symmetry of the initial lattice. The 
functions (7) are connected with the real basis, ipi oc d z i 
and ifi2 oc d x i_ y i : 

-01 = {<Pi +^ 2 )/V2 ; V2 = (Vi - «<P2)/\/2 (7') 

The tight-binding spectrum of the e2 S -electrons consists 
of two branches, which are discussed below. 

Before proceeding further, we write down explicitly the 
local, JT-part of the Hamiltonian (1) in the basis (7): 

~ 2 Qo ^ exp(-i(9) J (8) 

In the standard notations [7]: Q2 = Qo sm # ; Q 3 = 
Qo cos 9. Here Qo is the magnitude of the JT distor- 
tions, while the "angles", 9 — 0, 2ir/3, —2tt/2> correspond 
to elongations of the octahedron along the z, x, and y- 
axes, respectively. 

With the above in mind, consider the options for pos- 
sible ground states in manganites. In the case of stoichio- 
metric LaMn03 there is exactly one e2 9 -electron (Mn 3+ ) 
per unit cell. The ferromagnetic filling-up for two bands 
will produce a metallic spectrum. AF-ordering along one 
axis results in the same main contribution, — JhS (at 
large Jh), while the corresponding kinetic energy contri- 
bution becomes of the order of t 2 / Jh <C t. Second order 
terms in t in Eq. (4) imply that the spectrum is much 



less dispersive along the z-axis (M x = 0). Planes with 
the antiparallel core spins become almost isolated. 

AF-ordering along any other direction would lack elec- 
tronic transport along this direction, again shrinking 
bands. Thus, it is enough to have the AF order along 
one axis to access the major energy gain, — JhS. How- 
ever, to stabilize the AF-order along one axis, one needs 
to further reduce the total energy by the antiferroelastic 
JT-distortions along the remaining directions. 

Note again that after the AF-ordering sets in, say, 
along the z-axis, the communication between adjacent 
layers in real space becomes very weak. Therefore, at 
large Jh it is possible to reduce the problem, to main 
approximation, to the one of two-dimensional (2D) elec- 
trons. For decoupled layers, electrons inside the layer are 
ferromagnctically polarized along the 3 t2 9 -(core) spins. 
Again, in this 2D-problem there is one electron per (2D) 
cell. The metallic spectrum simplifies to the form: 

£ i^(p) = (A + B)(cosap x + cos ap v ) ± 
(A — B) { (cos ap x + cos ap y ) 2 + 

1 /2 

3(cos ap x — cos ap y ) 2 } (9) 

where 

i,k = (x,y,z) ; A oc ip ^z; x, y)tpi{z + a; x, y) , B oc 
if2(z; x, y)ip2{z + a; x, y) are the two overlap integrals be- 
tween the two sites separated by the lattice constant, a. 
Introduce now JT-distortions in the layer with the struc- 
ture vector (it, it) in a plane. If the JT-splitting is strong 
enough, then instead of (9), the new spectrum will consist 
of two pairs of bands separated by energy gaps over the 
entire (new) Brillouin zone, with the low-energy bands 
fully occupied. This completes the band picture for the 
insulating ^4-type ground state in LaMnC>3. 
For this interpretation one needs: 

Jh > t, gQ ; gQ > t (10) 

Eq.(10) docs not impose any severe limitations on the 
model parameters. The above physics should be present 
in any realistic band structure. Usually the band cal- 
culations use the experimental lattice parameters, and 
it is not straightforward to single out (e.g. in [3]) the 
competition of the effects considered above. 

Finding the single electron spectrum in presence of the 
JT-distortions from (8,9), writing down the total energy 
of the filled-up bands and subsequent minimization of the 
total energy, is a straightforward but rather tedious task. 
In the general case, there is no small parameter. We will 
consider numerical aspects elsewhere. 

Below we take advantage of the transparent physics 
of the independent 2D-layer approximation. First of all, 
the 2D-ferromagnetic ground state in each layer which 
was obtained in the mean field approximation, will be 
smeared out by fluctuations: ferromagnetism does not 
exist as a stable thcrmodynamical phase in 2D, un- 
less there is a weak interlayer coupling, t 2 j Jh, of eq. 



(6). This results in a rather low Neel temperature, 
T N ~ 150if, in LaMn0 3 , at f - O.leV, J ff - leV. 

To illustrate the above, consider the "symmetric" 
model [9], A = B = i/2 in (9): 

e(p) — t(cosap x + cosapy) = t(p) (11) 

Eq. (11) gives the two degenerate bands. For one polar- 
ized 2D-electron per unit cell, the two Fermi surfaces run 
along the lines p x ±p y = (±7r/a). In the presence of the 
JT-distortion (8) the new spectrum (in the reduced Bril- 
louin zone, with the structural vector (n/a, n/a)) consists 
of the four branches: £;(p) = ± || t(p) | ±A |, where 
A =| gQo/2 |. At A > 2t there are two pairs of energy 
bands, and the insulating gap sets in at p x = p y = 0, 
equal to A — 2t. The distortion gQ — g 2 /J e i is defined 
from the on-site problem. 

Turning back to the parameters {Ju,Jei,g,A,B) of 
our general model, connections between gQo,J e i and A 
and B are established by minimization of the total (elec- 
tronic and elastic) energies with respect to Q and M x 
(see eq. (5)). The deformations of octahedra {Qo) are 
known from the structural data (see e.g., in [3]). One 
may also expect that B <C A: tunneling between Mn- 
sites takes places through the shared oxygen ions, and 
the notion that the d-shell size, d, is small compared to 
the lattice constant, a(d <C a), is quite helpful, as oxides 
of transition metals are concerned [7, 13, 14]. Currently 
the spin wave spectrum in LaMn03 is available [15]. Al- 
though we postpone the study of the spin wave spectrum 
for the future, it is worth mentioning that coexistence of 
the localized 3 i2 S -spins and itinerant e2 S -electronic states 
(which in the AF-environment are not characterized by 
the spin-projection quantum number), may result in de- 
viations from an anisotropic Heisenberg model. 

Finally, we briefly discuss the concentration depen- 
dence. At small x the two mechanisms lead to local- 
ization of a doped hole. The first one is directly related 
to the 2D physics. According to (1), interaction of a hole 
with the JT-modes may decrease the hole's energy by an 
amount of order of g 2 /J e i- If its energy is now below 
the band's bottom, the hole becomes self-trapped. Two- 
dimensionality plays the crucial role in that there are no 
energy barriers for the process [12]. If a dimcnsionlcss 
parameter C <~ g 2 /J e iW (W is a band- width) exceeds 
a critical value of order of unity, no itinerant states are 
possible in the 2D system. The second mechanism is the 
Coulomb binding of holes near dopants [8] . 

At larger x the band approach becomes useless. In- 
stead, we adopt the point of view of the percolation the- 
ory. Let us start first with the cluster description [11]. 
In the so-called site problem, consider formation of an 
infinite cluster of the neighboring divalent ions. For the 
simple cubic lattice the critical x cr (s) = 0.31. This value 
is not universal. 

It is well known [11] that correlations between sites 
rapidly reduce x c (s) to the value x c ~ 0.16, which is also 
a threshold concentration for a continuous percolation 



problem. Although holes are located near divalent ions 
(forming the cluster's skeleton), the hole's hopping takes 
place along the adjacent Mn-sites correlated by the fact 
that the wave function of a single hole is spread over a few 
interatomic distances tending to suppress JT-distortions 
and align spins ferromagnetically in its vicinity. The clus- 
ter's structure is complicated, however, at x > x c one 
may expect that locally its properties are close to the 
ones of the homogeneous metallic ferromagnetic state. 
The rest of the sample does not conduct. 

Another plausible view is that the material may sep- 
arate into coexisting insulator and metallic phases [10]. 
In the cluster approach the Coulomb energy is playing 
the crucial role confining holes in the skeleton's vicinity. 
In the continuous description the Coulomb energy would 
limit sizes of the phases' domains. In both approaches the 
Coulomb energy may be strongly reduced by sufficiently 
large dielectric constant of the insulating phase and due 
to metallic screening in the ferromagnetic phase. 

With the temperature increase, concentration of the 
metallic phase decreases: temperature-induced JT- 
distortions [6] provide stronger localization of "doped" 
holes. From such a point of view the metal-insulator 
transition temperature, T* , may also be treated as a per- 
colation point (at fixed x), and this agrees well with the 
extreme sensitivity of CMR to applied magnetic field. 

To summarize, the view of LaMnC-3 as a band insula- 
tor, is consistent with the main experimental facts. Large 
intraatomic Hunds' exchange, Jh, forces the cubic sys- 
tem to order antifcrromagnetically along one direction 
because the band energy of electrons can then be effec- 
tively reduced by the JT-instability resulting in the new 
(tetragonal) lattice (y/2a x V2a x 2a). At Jh S> t 7 gQ n , 
the e2 9 -electrons form almost disconnected 2D-layers. 
This manifests itself in the low Neel temperature. The 
percolation analysis can explain the value of the critical 
concentration, x c ~ 0.16 and CMR itself suggesting new 
interesting physics. 
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